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Abstract 

Measurements of oscillation frequencies of the Sun and stars can provide important independent 
constraints on their internal structure and dynamics. Seismic models of these oscillations are used 
to connect structure and rotation of the star to its resonant frequencies, which are then compared 
with observations, the goal being that of minimizing the difference between the two. Even in the 
case of the Sun, for which structure models are highly tuned, observed frequencies show systematic 
deviations from modeled frequencies, a phenomenon referred to as the “surface term”. The dominant 
source of this systematic effect is thought to be vigorous near-surface convection, which is not well 
accounted for in both stellar modeling and mode-oscillation physics. Here we bring to bear the 
method of homogenization, applicable in the asymptotic limit of large wavelengths (in comparison 
to the correlation scale of convection), to characterize the effect of small-scale surface convection on 
resonant-mode frequencies in the Sun. We show that the full oscillation equations, in the presence of 
temporally stationary 3-D flows, can be reduced to an effective “quiet-Sun” wave equation with altered 
sound speed, Briint-Vaisala frequency and Lamb frequency. We derive the modified equation and 
relations for the appropriate averaging of three dimensional flows and thermal quantities to obtain the 
properties of this effective medium. Using flows obtained from three dimensional numerical simulations 
of near-surface convection, we quantify their effect on solar oscillation frequencies, and find that they 
are shifted systematically and substantially. We argue therefore that consistent interpretations of 
resonant frequencies must include modifications to the wave equation that effectively capture the 
impact of vigorous hydrodynamic convection. 

1. INTRODUCTION 

Measurements of oscillation frequencies of stars pro¬ 
vide some of the strongest constrain ts on stellar structure 
models. Rece nt missions Kepler (iBorucki et al.l I20IOI1 
and COROT (|Auver6fne et alJ 12009( 1 have provided us 
with high quality asteroseismic data, which have revealed 
structural features of Sun-like stars in unprecedented 
detail. There have been several successful solar mod¬ 
els which have reproduced the oscillation frequencies to 
within a few parts in thousand. Despite this degree of 
accuracy, the modeled frequencies display systematic de¬ 
viations from the observed ones. Understanding the root 
of these deviations is important both from theoretical 
and observational point of view — not only is it a win¬ 
dow into the missing physics, it is also a hindrance in 
htting stellar models. With the ever increasing precision 
in measured frequencies, it is imperative to understand 
the physical reason behind the trend and to correct for it 
while extracting information about the Sun at the same 
time. 

these shifts are ob- 
model and modal ef- 
Some of the possible 

sources for this bias are: 

1. Errors in physics that go into the modeling, includ¬ 
ing turbulent-pressure contributions (model effect). 

2. Non-adiabaticity prevalent in the outer layers of 
the Sun (model+modal effect). 

3. Imprecise modeling of propagation of waves 


There are various reasons why 
served, which were categorized as 
fects by ([Rosenthal et al’l 119991) . 


through convective flows near the surface of the 
Sun (modal effect). 

4. Non-linearity in the wave equation (modal effect). 

It is unlikely that any one of these mechanisms would ex¬ 
plain the observed deviations completely, rather a com¬ 
bination of these corrections would lead to a better es¬ 
timate of oscillation frequencies, and possibly eliminate 
the observed bias. We would like to differentiate between 
the first and third points — the first assumes a hydro- 
dynamic equilibrium model, but treats oscillations as if 
they were about an equivalent hydrostatic background, 
while the third treats the oscillations about a dynamic 
background. In this article we focus on the third as¬ 
pect, and propose a method to correct the modal effect. 
Strong convective flows with speeds close to the sound 
speed exist near the surface of the Sun. We show that 
advection of traveling waves by these flows can have a sig¬ 
nificant impact on the resonant mode frequencies, and it 
is therefore important to account for this effect in order 
to extract the correct physics from measurements. 

Acoustic modes {p modes) are standing waves bounded 
from below by a lower turning point set by the Lamb 
frequency, and from above by an upper turning point 
set by the acoustic-cutoff frequency. The upper turning 
point for modes with frequencies below 2 mHz is signif¬ 
icantly below the surface of the Sun, but modes with 
frequencies lying in the range 2 mHz to 5 mHz encounter 
an upper turning point close to the surface. These are 
the modes for which the measured frequencies display a 
systematic deviation from model, hence it seems likely 
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that the missing physics lies close to the surface, and the 
deviation is often termed ‘surface effect’. The deviation 
is predominantly a function of the frequency of a mode, 
hence corrective terms should also be functions of the 
mode frequency. An early attempt at m itigating the fre¬ 
quency differences was by iGoughI (jlQQOtl , whose prescrip¬ 
tion was in the form of a power law correction based on 
the physical process involved: fibrillated magnetic fields 
lead to a change in sound speed, and hence an observed 
frequency v would be shifted hj Sv (x where I 

is the mode inertia; while alternate para metrizations of 
convec tion lead to a shift Si/ oc 1/ {vT). iKieldsen et al.l 
(j2008ll had studied the shifts in the Sun, and with an 

empirical power law model Sv = aiyjv^f ^ were able to 
correct f or the shifts in a Ce n A, a Cen B and /3 Hyi. 
Recently iBall fc GizonI (120141 1 have shown that a combi¬ 
nation of the two effects studied by Gough — the cubic 
shift and the inverse shift — is able to fit the deviations 
better than Kjeldsen’s power law. These corrections are 
satisfactory from a model-fitting viewpoint, but by elim¬ 
inating the frequency differences altogether we lose in¬ 
formation about their origin. 

Despite the high flow speeds, there still exists a small 
parameter in the wave equation — the granulation length 
scales are usually much smaller than horizontal wave¬ 
lengths of lo w-£ modes, whi c h are of the order of the 
solar radius. iHanasoge et al.l (j2ni.‘l[ l had designed an al¬ 
ternate approach to deal with small-scale sound speed 
perturbation, by expanding the wave equation in terms 
of the ratio of the two length scales, followed by averag¬ 
ing over small scales. They had shown that if the hori¬ 
zontal wavelength of a propagating wave is much larger 
than the length scale over which the background medium 
fluctuates, it is possible to replace the background by a 
homogeneous, coarse-grained one. In this paper, we use 
a similar analysis in presence of small-scale flows, and 
derive an expression for an effective speed. A way to 
think of homogenization is as follows: near-surface con¬ 
vection is a stochastic process, evolving both temporally 
and spatially. Here we assume temporally stationarity, 
with waves propagating through a frozen configuration of 
convective flows. Precise measurements of resonant fre¬ 
quencies require taking spectra of long sequences of raw 
observations of the wavefield — at a minimum, several os¬ 
cillation periods but typically much longer than that, on 
the order of a million periods, so as to resolve the mode 
line-width and beat down the noise sufficiently. Over this 
time, waves will have propagated through a large num¬ 
ber of realizations of convective flows. Thus one may 
think of this as an ensemble of wave equations, each as¬ 
cribed a set of coefficients corresponding to different re¬ 
alizations of the process of convection, drawn from the 
same probability density function. Each realization of 
convective flows corresponds to a wave solution, and the 
question then becomes: how to obtain an ensemble av¬ 
erage of these wave solutions? Homogenization provides 
an answer to this question, specifically in the limit where 
the horizontal wavelength is much larger than the length 
scale corresponding to the flows. The principle behind 
homogenization is that in cases where waves and flows 
are scale-separated, measurables such as frequencies can 
be reproduced by studying wave-propagation through a 
quiet, homogeneous medium obtained by appropriately 


averaging over small-scale flows. The method has an ad¬ 
vantage that the effective speed which governs wave be¬ 
haviour on a large-scale emerges naturally out of the av¬ 
eraging process. In the asymptotic limit, this is the speed 
which shall appear in the dispersion relation, and hence 
shall determine the observed frequency differences. The 
method of homogenization allows an extension to ran¬ 
dom media where the distribu tion function is translation- 
ally i nvariant and bounded (I Papanicolaou &: VaradhanI 
Il982f) , and thus holds promise in progressing beyond our 
simplified analysis to real solar models. _ 

Several au t hors s uch as iMurawski fc Roberts! (jI993f ): 
iDuvall et al.l (|I998f) had studied the effect of random 
scatterings by convective flows on surface gravity waves 
(/ modes), arriving at corrections to the dispersion rela¬ 
tion in order to match the observed frequencies and line- 
widths. A different treatment to captur e the effect o f 
flows on p modes had been carried out bv iBrownl (jI984[l . 
who had proposed a correction of the form (l — (uf)/c^) 
by carrying out a horizontal average over the wave equa¬ 
tion. This form of the correction comes closest to our ap¬ 
proach, which generalizes the analysis to arbitrary small- 
scale vector flow velocities. An asymptotic analysis of p- 
mode frequen cy shifts for small Mach n umbers had been 
carried out bv iStix fc Zhugzhdal (1200411 , but whether the 
analysis extends to transonic flows near the solar surface 
remains unclear. 

2. GOVERNING DYNAMIGS 

Waves on the Sun are modeled as small-amplitude per¬ 
turbations about a background, which we assume to be 
in a steady state. The background medium is charac¬ 
terized by its pressure p (x) , density p (x) , flow velocity 
u(x), and gravitational potential (/)(x). Waves shall in¬ 
duce oscillations about the temporarily stationary back¬ 
ground, and the velocity of a moving blob would be given 
by a sum of the background flow velocity and the mode 
velocity. We work with the vector flow velocity u (x) 
without splitting it up into horizontal and vertical com¬ 
ponents, since both the components repeat over small 
scales horizontally and hence can be treated analogously. 
In further analysis, we shall suppress explicit spatial and 
temporal dependence for notational convenience. The 
full magnetohydrodynamic equation for wa ves traveling 
throug h non-uniform flows was analyzed by iWebb et al.l 
(l2f)n^ . we replicate the equation of motion here. The 
background flow field u satisfies the equation 

V • (puu-l-pl) -I- pV(/) = 0, (I) 

where I denotes the second-rank identity tensor dehned 
as lij = Sij. In presence of a flow field u, the material 
derivative of a small-amplitude displacement ^ (x, t) is 
defined as 

i = U + u- VC (2) 

We assume that there exists a flow field which satishes 
the background equations, and concentrate on the wave 
equation. The acoustic wave equation, to linear order in 
the perturbations, is given by 

dt (pi) + V • {pui + I{p- pc^) V-^-p {VCf } 

+pi • VV</) = 0, (3) 










































3 


where T in the exp onent denotes the transpose operator 
(|Webb et alJl2005r) . We shall stick to the Cowling ap¬ 
proximation as our focus is on correcting for the wave 
speed; the analysis can be extended to a general case in 
an analogous manner. 

A propagating wave is scattered out in all direction by 
the background flow field, which results in a change in 
the phase of the forward-traveling wave. Multiple scat¬ 
tering would lead to repeated phase shifts, leading to an 
altered phase velocity of the wave. This modified ve¬ 
locity would govern the dispersion relation, which would 
in turn determine the oscillation frequencies. In order 
to model this scattering perturbatively, we resort to the 
multiple scale approach of spatial homogenization. 

3. SPATIAL HOMOGENIZATION 

Low-degree resonant modes on the sun have horizontal 
wavelengths (A) comparable to the solar radius, which 
is much larger than the typical length scales of granula¬ 
tion on the Sun (L). In such a case, the dynamics can 
be parametrized in terms of the scale ratio e = L/X. 
In the limit e -C 1, we can replace the rapidly fluctuat¬ 
ing medium by a homogeneous, albeit anisotropic one. 
To illustrate this approach, we study a special case of a 
periodic flow in a Cartesian box, where we denote the 
vertical depth by the coordinate z, and label a point on 
a horizontal plane by the coordinates {x,y). The back¬ 
ground medium repeats after a length L along both the 
horizontal directions x and y, and is stratified in the ver¬ 
tical direction. On a particular horizontal plane, the flow 
takes the form of tessellated square cells. To study the 
dynamics in the asymptotic limit £ <C 1, we introduce 
two spatial coordinates which are scale-separated: a slow 
coordinate x = {x, y, z) and a horizontal fast coordinate 
x'^ = {xje,yle). Note that the vertical component of 
mode wavelength can be comparable to the scale of flows, 
so the advection terms might not be scale-separable ver¬ 
tically. Bearing this in mind, scale separation is being 
carried out strictly in the horizontal directions, labeled 
by a; — y or equivalently by xj^. The depth coordinate 
z is common to the flows and the wave, and it charac¬ 
terizes the dynamics completely in the vertical direction. 
In terms of these new coordinates, the spatial gradient 
can be rewritten as 

V = Vx -f I/eVh, (4) 

where Vx acts over the coordinates denoted by x, and 
Vh acts over the fast horizontal coordinates denoted by 

<• 

We expand the wave field ^ as a perturbative series in 
the small parameter e as 

^ (x, x'^, t) = lo (x, X^,, t) -f £|i (x, x;,, t) 

+e ^^2 (x, x'^, t) + 0 (e^) . (5) 

The term at order e° represents the large-scale charac¬ 
teristics of the wave field, while the higher-order terms 
arise as the incoming wave is scattered by the flows. This 
approach differs from standard perturbative expansions 
in two ways — firstly the expansion parameter is not in 
terms of the flow speed u but the length scale L; secondly 
the zeroth-order field is not a solution to an equation 
without flows, but senses an average effect of the flows. 
Our approach would be to replace the spatially varying 


medium by an equivalent quiet background, which is dif¬ 
ferent from what would have been there had there been 
no flows. We substitute Equation (O in Equation ([3)) 
and use Equation ([2]) to write down the equations order 
by order in e. 

At order we get 

+ pIo • VhVh0 = 0, (6) 

where 

£(-2)| = Vh-{puu.Vhl+I {p - pc^) Vh-I-P y 

We introduce an alternate definition of the operator 
£(“ 2 ) ] 3 y defining the transpose map explicitly. We in¬ 
troduce the following tensors and tensorial relations: 

• The double-dot product of two dyadic tensors A 
and B defined as A : B = AijBij, where repeated 
indices are summed over. This definition can be 
extended to higher order tensors A and B, where 
the last two indices of A and the first two of B 
shall be contracted. 

• A fourth order tensor T defined as Tyfcz = SuSjk- 
The tensor T has the property T : D = for a 
dyadic tensor D. The definition can be extended to 
higher-order tensors, where the transpose operator 
T flips the first two indices. 

• A fourth-rank tensor I4, defined as 14 ^^; = 6ikSji^ 
that satisfies I4 : D = D for any dyadic tensor D. 

We define the tensor 

f = puu • 14 -|- I (p — pc^) I — pT ( 8 ) 

in terms of which we can rewrite the operator as 

£(-2)| = Vh-(f :Vh|). (9) 

We shall assume henceforth that the gravitational poten¬ 
tial (j) is independent of the fast coordinate x^, which lets 
us drop any term like or VhVxf/*- This assumption 
is reasonably justified, since density in the near-surface 
layers is much smaller than the interior and hence the 
gravitational potential would depend only weakly on the 
flows near the surface. Under this assumption, Eq. 
becomes 

= 0. (10) 

A particular solution to Equation ITOl) is a zeroth-order 
field ^0 that is independent of the fast coordinate xj^. 

This is a sufficient, but not a necessary condition. A 

general solution to Equation (fTOl) can be a function of xj^, 
which means the operator might be singular. We 

shall need to invert this operator at order , so it might 
be necessary to define the inverse precisely by imposing 
boundary conditions on the fields, or as a generalized 
inverse, which is beyond the scope of this paper. We shall 
proceed assuming that ^0 is independent of xj^, which fits 
into our picture of perturbations about a homogeneous 
background. 

Collecting terms of the order e“^, we obtain 

-f = 0, 


( 11 ) 
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where 

£(-i)| = Vh- (f : Vxl), (12) 

and is as defined in Equation ([HI). We use the fact 

that is independent of to rewrite Equation m as 

Vh • (f : Vhli) = - (Vh • f) : Vx (13) 

We note two facts - firstly, Equation (fT51) is linear in 
^ 1 , secondly, the right hand side is a product of two 
terms, one of which varies horizontally over the slow 
scale X while the other varies over the fast scale xj^. 
This motivates us to seek solutions for by separat¬ 
ing the slowly varying horizontal components of x from 
the rapidly varying components of x)j. We substitute 

|i=J(x;,z): Vxlo (14) 


resembles a wave speed squared term that governs wave 
propagation on large s cales. This is similar in form to 
the result obtained by iBrownI (Il984f) . except the term 
J which arises from the separation of scales. We can 
substitute for F from Equation ((Sj) to recast Equation 
(1^^ in a more tangible form as 


C = 




II-uu L-f-T 

P 


-\- terms arising from cell problem. 


(23) 


Convective cells on the surface of the Sun have a struc¬ 
tural form where hot gas rises up near the center and 
cooler gas sinks down around the edges of the cell. We 
assume that the flow in a particular cell is anti-symmetric 
about the cell center, hence (pu) = 0. This lets us drop 
the advective term in Equation (1211) and obtain 


where J (xj^, z) is a third-rank tensor. Substituting Equa¬ 
tion (fT^ in Equation (fTT|) . and using the fact that is 
independent of x'^, we obtain 


£(-2) J + Vh • f 


: Vx ^0 = 0 . 


(15) 


We seek a solution J which satisfies 

= -Vh • f. (16) 

This is conventionally known as a cell problem (see e.g. 
iBensoussan et aHllfiTSf ). 

At order e° we get the equation of motion for the 
zeroth-order field 

pdtio + pn-dt (Vxlo + Vhli) 

+p|o-VxVx<)) = 0, (17) 

where 

/:(o)| = Vx-{f : Vxl}, (18) 

/:''’l = Vx-{f : Vhl}. (19) 

We know that the first order field depends on the 
zeroth-order field through Equation (IT^ : we substi¬ 
tute this in Equation (flT)) and average over the fast co¬ 
ordinate x(j. The averaging process restricts us to the 
zero-spatial-frequency mode, however we do retain infor¬ 
mation about smaller length scales through J. Given 
a horizontal cell B — [0, L) x [0 ,L), we define the cell 
average of a function / (xj^) as 


(p)92|o-Vx- 


(pC) : Vxlo 


+lo-(pVxVx(?!>) = 0 . 


(24) 


Equation (I24p is the homogenized wave equation, 
wherein we have replaced the spatially fluctuating back¬ 
ground medium by an effective homogeneous, anisotropic 
one. The effective speed depends on both the horizon¬ 
tal and vertical flow profiles, the sound speed profile and 
the density. If density variations in a cell are spatially 
correlated with the flow speeds, a description in terms of 
the average speeds and average densities which ignores 
correlations might not be sufficient. 

Periodicity of flows along specific directions leads to 
a tensorial wave speed, hence we expect anisotropy in 
wave propagation. Equation (1^^ is a standard eigen¬ 
value problem which we can solve to find the dispersion 
relation w (n, k) for vertical order n and horizontal wave 
vector k. This dispersion relation captures the extent 
to which the frequency of a monochromatic wave would 
appear to differ from the canonical dispersion relation 
Lo = ck, which is governed only by the adiabatic sound 
speed. As an example, we consider waves propagating 
in a 2-dimensional uniformly dense fluid. We assume 
that the fluid film is tiled with small-scale, periodic flows 
u = (ux, Uy). A plane wave propagating along the x axis 
would have the form 


lo (x) = Aexp(i(fc:j,a;-a;f)), (25) 

where A = (Ax^Ay) is the amplitude, kx is the x com¬ 
ponent of the wave-vector and ui is the frequency of the 
wave. Neglecting changes in gravitational field 'Vx4> in 
the plane, we can rewrite Equation (I24|) in Fourier space 
as 


= ( 20 ) 

The terms and in Equation (TTfl) are to¬ 

tal divergences in xj^, so they drop off on averaging as 
we impose periodic boundary conditions. We obtain 


{p)dt^o + 2(pu) • Vx dt |o 


-Vx 


{pC) : Vxlo 


+ lo ■ (pVxVxt))) 


0 , ( 21 ) 


where 

C = f : (14 + Vhj) ( 22 ) 


Ax 

A,, 


= ki 


(c^ -ul+ C.pi) (C.P 2 ) 

(C.ps) {-ul + C.p4) 


Ax 

Ay 
(26) 

where the terms referred to as c.pi to c.p 4 arise from the 
solution to the cell problem. The corresponding disper¬ 
sion relation can be obtained by diagonalizing the matrix 
in the second term of Equation (I^Hl) . If the terms c.p^ to 
c.p 4 are negligible, the dispersion relation would have the 
form u) = (y^c^ — Ux)kx, but this might not be a good 
approximation in general. 

We note that Equation (IMl) is not damped at the lead¬ 
ing order, so attenuation, if any, would be captured at a 
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higher order. This is consistent with observations, since 
the observed p-mode line-widths are much smaller than 
the corresponding frequencies. 

4. NUMERICAL TEST OF FLOW HOMOGENIZATION 

In this section we take a simple example of waves prop¬ 
agating through flows in a 2-D plane, and show that we 
can effectively reproduce the wave speed by homogeniz¬ 
ing over small scale flows. A thorough analysis of the 
approach would include realistic pressure, density and 
sound-speed profiles, possibly starting from simulations 
of near-surface structure of stars. In this work, we fo¬ 
cus only on sound-speed variations induced by flows, 
and assume that the background pressure and gravita¬ 
tional field are unchanged in presence of the flow. We 
also assume that the background medium is uniformly 
dense throughout. Under these assumptions. Equation 
m simplifies considerably as we can drop the terms 

V- ^IpV ■ ^ — p and Wc^. We combine Equa¬ 

tions @ and ([3]) to obtain 


as yet arbitrary, wave-vector, and the function rjk has 
the periodicity of the flow. We study waves propagat¬ 
ing along the x-axis which have ky = 0. We impose 
a large-scale periodicity condition | (cc-I-27ri?0, j/; w) = 
^{x,y;uj), which restricts us to discrete values of the 
wave-vector kx = {27 :/Rq) n for integral n. Correspond¬ 
ing to each value of n we obtain an eigenvalue problem 
from Equation (1301) . which we solve numerically for the 
frequency w {kx)- We refer to this relation as ‘true dis¬ 
persion relation’. In the limit of infinite wavelength, we 
can replace the fluctuating medium by a homogeneous 
one, given by Equation (1331) . We make one further ap¬ 
proximation — because of potential non-uniqueness in 
computing J, we restrict ourselves to the zeroth-order 
homogenized speed given by 


C = c2ll-uu-i4. (31) 


The homogenized equivalent to Equation (1331) in this case 
is 


- Vx • 


(C) : Vxlo 


= 0, 


(32) 


df$, + 2u • Vdt$. + (u • V)^ I - V (c^ V • I) = 0. (27) 

We assume that the waves are propagating in a 
medium that is periodic along both the directions; in 
other words, the unit box B — [0, L) x [0, L) is tessellated 
to hll up the entire two-dimensional plane. We choose to 
work in Cartesian coordinates {x,y), with the origin at 
the center of a cell. As an example of a flow-field in a 
cell, we choose a vortex centered about the cell-center. 
The functional form of the flow is given by 


{xy - yx) 

u = uq -exp 


{x^ + y^) \ 
2 ct 2 j ’ 


(28) 


where uq is a velocity scale, cr is a length scale of the 
flow related to the cell-size L. The sound speed c would 
depend on several factors such as the pressure in a cell 
and the equation of state. For our illustrative purpose, 
we choose a simple model given by 

c {r)'^ =clo + {4 - c^) exp ^ , (29) 


where C is as defined in Equation (|3T|) . We compute the 
dispersion relation corresponding to Equation p2[) . and 
refer to it as ‘homogenized relation’. It is worth pointing 
out that the approximation made in Equation m is 
crude; one can not neglect J if Mach numbers are close 
to 1, and it might be necessary to compute J by solving 
Equation (1331) to achieve a higher level of accuracy. 

We compare the zeroth-order homogenized dispersion 
relation with the true dispersion relation in Figure [TJ 
We also compare it with dispersion relations of the form 
uj = {c)k and lo = y/JcFjk, which do not explicitly ac¬ 
count for flows. We have listed frequency differences for 
several different parameter values in Table [31 We find 
that the homogenized relation replicates the true disper¬ 
sion relation better than the alternate relations. The rel¬ 
ative frequency differences corresponding to Run 2 have 
been plotted in Figure [3] Comparing with Run 1, we 
note that the frequencies predicted by the zeroth-order 
approximation to the homogenized dispersion reduce in 
accuracy as the Mach number increases. 

5. IMPACT ON SOLAR FREQUENCIES 


where cq is the sound speed at the center of the cell and 
Coo is the sound speed far away from the center of the 
cell. We choose a cell length-scale L — 1 Mm, and values 
of sound speed and flow parameters as listed in Table [H 
where u refers to the magnitude of the flow velocity u. 

We recast Equation (|37|) in the form of an eigenvalue 
problem by using Equation ([3]) along with Equation ([371) . 
Assuming harmonic time dependence, we rewrite the sys¬ 
tem of equations as 


— LUJ 



I(-u-V) I 
D I(-u V) 



where the operator D encapsulates the restoring force on 
the wave field ^ and is given by = V (c^V • ^). We 
note that both u • V and D are periodic with the same 
length scale as the flow, hence we use Bloch’s theorem 
to reduce the problem to one cell. According to Bloch’s 
theorem, the wave-field ^ can be written as ^(x;a;) = 
exp (ik • x) rjk (x;a;), where k = {kx, ky) is a horizontal, 


Correcting the solar oscillation equations involves cor- 
re cting for both modal a nd model effects, a s described 
bviRosenthal et aP (jl999ll . In a recent paper, iPiau et id] 
(|20I4f) have explored the impact of various parametriza- 
tion of stellar convection — from mixing-length theory 
([Spruitl 1197711 to patched hydrodynamic models — on 
oscillation frequencies. They have used horizontally- 
averaged sound speeds from these models in the wave 
equation. We have shown that the presence of flows in¬ 
troduces further advective corrections in the wave speed, 
and the horizontally-averaged sound speeds do not rep¬ 
resent the true wave speed. The actual speed of the 
waves is obtained from Equation (1241) , and the difference 
between this speed and the horizontally-averaged sound 
speed can be substantial near the solar surface. 

In this section, we start from Model S 
([Christensen-Dalsgaard et al.l 1199611 which uses the 
mixing-length formulation of convective flux. We correct 
wave speeds using flow profiles obtained from simu¬ 
lations by Schiissler (personal communication) using 
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(jj = ^{c^)k. We see that the homogenized relation reproduces mode frequencies with higher accuracy than the other relations in the low 
frequency limit. 
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Fig. 2.— Relative frequency differences as obtained in Run 2, where the markers represent the same quantity as Figure [T] Run 2 
corresponds to a higher Mach number than Run 1, and we find that the 0-th order homogenized approximation is not as effective at 
reproducing the true dispersion relation as in Run 1. 
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CO 

Coo 

(c) 

'^^max 

(7 

(u) 

Run 1 

10 km/s 

11.7km/s 

11.6 km/s 

4 km/s 

0.1 Mm 

0.5 km/s 

Run 2 

10 km/s 

11.7km/s 

11.3 km/s 

4 km/s 

0.2 Mm 

2.0 km/s 

Run 3 

10 km/s 

13.7km/s 

12.8 km/s 

6.2 km/s 

0.2 Mm 

3.0 km/s 

Run 4 

10 km/s 

11.7km/s 

10.7 km/s 

4 km/s 

0.4 Mm 

3.5 km/s 


TABLE 1 

Sound speed and flow speed corresponding to two runs with different flow speeds in a cell. 



(Run 1) 

(Run 2) 

(^'^)avE (Run 3) 

(Run 4) 

to — Chom^ 

0.009 /^Hz 

0.07/iHz 

—0.5 ^Hz 

0.29 ^Hz 

iO = {c)k 

—0.15 fiHz 

—0.59 ^Hz 

— 1.8 ^Hz 

— 1.36 fiHz 


—0.16 fiHz 

—0.62 ^Hz 

— 1.9 ^Hz 

— 1.40 fiHz 


TABLE 2 

Comparison between the frequency differences corresponding to the different dispersion relation, averaged over the 
range 4/^Hz <u< 120/iHz. The homogenized relation, which accounts for flows explicitly, comes closest to the true 

DISPERSION RELATION. 





the MURaM code (|Vogler et alJ 120051) . The flow-speed 
profiles have been plotted in Figure [S] We start from 
equation (|2^ and work with low-degree modes (£ < 3), 
which are predominantly radial. We assume that the 
Mach number is low enough for the zeroth-order ho¬ 
mogenization approximation (Equation I31|l to hold, and 
retain only the dominant component of C which goes 
roughly as (? — . The exact form of u which goes into 

the correction depends on the direction of advection at 
every depth; for our purpose we shall restrict ourselves to 
the special cases of radial flow (u = uy) and tangential 
flow (u = ux). We refer to this modified model as Model 
Horn, and compute the eigenfrequencies corresponding 
to this model assuming the same boundary conditions as 
Model S. We compare these frequencies with observed 
solar frequencies obtained b y the Birmingham Solar 
Oscillation Network (BiSON) (|ChaDlin et aLll^OOTf) . We 
find that the presence of flow advection in the wave 
equation leads to a reduction in mode frequencies, which 
reduces the observed difference between observed and 
Model S frequencies. 

To study the impact of modal versus model correc¬ 
tions on frequencies, we compare the M odel Horn fre¬ 
quen cies with another hydrostatic model (jBasu fc Antial 
11994 henceforth called Model BA), which uses 
the Canuto-Mazzitelli form ulation of convective flux 
(|Canuto fc Mazzitellil II 99 HI . The frequency differences 
between these models and the observations by BiSON 
have been plotted in Figure IH We find that the modal 
frequency shifts due to the presence of advective flows is 
comparable to those obtained from modeling differences. 
Both the effects need to be looked into while fitting solar 
frequencies. 

Some of the assumptions made in this section, such as 
diagonally-dominated and isotropic C, are expected to 
be violated considerably near the surface, but the present 
example demonstrates that the effect of flow correction 
can be substantial. A detailed computation, starting ei¬ 
ther from a fully hydrodynamic background or a patched 
hydrodynamic background, is necessary to make quanti¬ 
tative statements about frequency shifts due to convec¬ 
tive flows. 


6. CONCLUSION 

Observed solar oscillation frequencies can differ from 
modeled frequencies because of an inaccurate back¬ 
ground model, as well as incomplete description of mode 
physics. The presence of near-surface convective flows 


on the Sun results in mode frequencies which are dif¬ 
ferent from ones computed starting from a quiet back¬ 
ground. Three-dimensional convection simulations gen¬ 
erally result in an increased size of the acoustic cavity , 
thereby reducing the frequencies (jRosenthal et al.l 119991 
iPiaii et nil2fll4l l. In this paper, we have shown that the 
impact of advection on waves might result in frequency 
shifts of similar magnitudes. In presence of flows, the 
appropriate wave speed is not the horizontally-averaged 
sound speed, rather it is a combination of sound speed 
and an appropriate projection of the flow speed in the 
direction of propagation. 

Granulation length scales on the solar surface are 
much smaller than horizontal wavelengths correspond¬ 
ing to low-£ modes, which lets us replace the irregular 
background medium by an effective, homogeneous one. 
Acoustic waves propagating through steady background 
flows experience a change in phase speed on being scat¬ 
tered by the flows. The effective wave equation in the 
large-wavelength approximation has the form 


{p)d^^o - V • 


(pC) : V|o 


+ lo ■ (pWcf)) = 0, 


where the effective speed C is a tensor, given by Equation 

( 1221 ). 

The impact of time evolution of flows, missing in the 
present analysis, needs to be studied in detail. The spa¬ 
tial scales of waves and background flows are signifi¬ 
cantly different, thereby allowing a two scale analysis. 
The timescales of granule evolution are, however, not 
very different from time periods of the wave. Overlap of 
scales generally results in strong scattering and forcing, 
thus a picture of flows evolving in time might be nec¬ 
essary to capture frequency shifts and damping. Added 
to this are the impact of super-adiabaticity and mag¬ 
netic fields, which have been ignored in the present anal¬ 
ysis, but might be important in governing the spectrum. 
Nevertheless the present analysis provides a useful start¬ 
ing point to model wave propagation through convective 
flows and resultant frequency shifts. 
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Fig. 3. — Root-mean-square flow speeds obtained from three-dimensional simulation of the near-surface layers of the Sun flVogler et al.l 
120051) . 



BiSON frequencies Hz) 


Fig. 4.— Eigenfrequencies corresponding to Model S (IChristensen-Dalsgaard et al.lfT996l . circles) compared with the modified models. 
The squares indicate a model where advection is predominantly radial and the st ars indicate a model where advection is predominantly 
tangential. The dotted line represents eigenfrequencies corresponding to Model BA JBasu fc: Antialll994l) . We find that the modified models 
display frequency shifts comparable to Model BA, which shows that modal effects can be significant near the surface and need to be taken 
into account. 
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